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Abstract 

We present a full next-to-leading order code to calculate the photoproduction of prompt 
photons. The code is a general purpose program of "partonic event generator" type with large 
flexibility. We study the possibility to constrain the photon structure functions and comment 
on isolation issues. A comparison to ZEUS data is also shown. 
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1 Introduction 



High energy electron-proton scattering at the DESY ep collider HERA is dominated by photopro- 
duction processes, where the electron is scattered at small angles, emitting a quasireal photon which 
scatters with the proton. These processes are of special interest since they are sensitive to both 
the partonic structure of the photon as well as of the proton. In particular, they offer the possibil- 
ity to constrain the (presently poorly known) gluon distributions in the photon, since in a certain 
kinematical region the subprocess qg — > 75, where the gluon is stemming from a resolved photon, 
is dominating. Up to now, the experimental errors were too large to discriminate clearly between 
different sets of gluon distributions in the photon, but a high statistics analysis of the 1996-2000 
HERA data on prompt photon photoproduction announced by the ZEUS collaboration will shed 
new light on this issue. 

The calculation of higher order corrections to the Compton process 7(7 — ► jq has been initiated 
some time ago |f§> The most recent calculations for prompt photon photoproduction have 
been done by Gordon/ Vogelsang || for isolated prompt photon production, Gordon for photon 
plus jet production and by the group Krawczyk/Zembrzuski Q for both the inclusive case and 
7+jet. However, all of these calculations contain certain drawbacks. In isolation is implemented 
by adding a subtraction term evaluated in the collinear approximation to the fully inclusive cross 
section. The programs of and Q do not contain the full set of NLO corrections. In R], those 
parts where the final state photon comes from fragmentation of a hard parton were included only at 
leading order, arguing that isolation cuts will suppress the fragmentation component in any case to a 
large extent. Moreover, the box contribution has not been included. In ]8), higher order corrections 
are included only for the case where initial and final state photons are both direct. So not only the 
contributions from fragmentation, but also the case where the initial photon is resolved are included 
at Born level only. However, the box contribution has been taken into account. 

The calculation presented in this paper takes into account the full NLO corrections to all four 
subparts. The corresponding matrix elements already have been calculated and tested in previous 
works ||, [l(|. A major advantage of the present code is also given by the fact that it is constructed 
as a "partonic event generator" and as such is very flexible. Various sorts of observables matching 
a particular experimental analysis can be defined and histogrammed for an event sample generated 
once and for all. This strategy already has been applied to construct NLO codes for 77 production 
(DIPHOX) and one or two jets photoproduction [ fl2"| . 

The paper is organized as follows. In section 2 we first describe the theoretical framework and 
the treatment of the infrared singularities. Then we discuss the implementation of isolation cuts 
and outline the structure of the code. Section 3 is devoted to phenomenology. We study the effect 
of isolation, determine the kinematic region which is most sensitive to the gluon distribution in the 
photon and illustrate the sensitivity of the cross section to the energy of the incoming photon. We 
give results for inclusive isolated prompt photon production and compare with a recent analysis of 
ZEUS data before we come to the conclusions in section 4. 

2 Theoretical formalism and description of the method 

In this section the general framework for prompt photon photoproduction will be outlined. We will 
review the contributing subprocesses, the treatment of infrared singularities and the implementation 
of isolation cuts. 

2.1 The subprocesses contributing at NLO 

The inclusive cross section for ep — > 7 A" can symbolically be written as a convolution of the parton 
densities of the incident particles (resp. fragmentation function for an outgoing parton fragmenting 
into a photon) with the partonic cross section a 
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da ep -'^ x {P p ,P e ,P J ) = Yl J d *ef dx pJ dzF a/e (x e ,M)F b/p (x p ,M p ) 

a.b.c 

da ab ^ cX (x p P p , x e P e ,P 7 /z, fi, M, M p , M F ) D j/c (z, M F ) (1) 

where M, M p are the initial state factorization scales, Mp the final state factorization scale and /j, 
the renormalization scale. 

The subprocesses contributing to the partonic reaction ab — > cX can be divided into four cat- 
egories which will be denoted by 1. direct direct 2. direct fragmentation 3. resolved direct 
4. resolved fragmentation. The cases "direct direct" and "resolved direct" correspond to c = 7 and 
Z? 7 / C (z, Mp) = <5 C7 (5(1 — z) in (|l|), that is, the prompt^] photon is produced directly in the hard 
subprocess. The cases with "direct" attributed to the initial state photon correspond to a = j, with 
Fry/ e approximated by the Weizsacker- Williams formula for the spectrum of the quasireal photons 

mv) = ^ll±il^ ln 9k4^1>^J^ll\. (2) 

7 2vr [ y miy 2 y J 

The "resolved" contributions are characterized by a resolved photon in the initial state where a 
parton stemming from the photon instead of the photon itself participates in the hard subprocess. 
In these cases, F a / e (x e , M) is given by a convolution of the Weizsacker- Williams spectrum with the 
parton distributions in the photon: 

F a / e (x e ,M)= dydx 1 f^(y)F a/l (x 1 ,M)8{x 1 y - x e ) (3) 
Jo 

Examples of diagrams contributing at Born level to the four categories above are shown in Figs. [[] 
and ||. 

In the case of the "direct direct" part, only the Compton process 79 — > jq contributes at leading 
order, at NLO the 0(a s ) corrections from 75 — > ^yqg resp. 75 — > ^qq and the corresponding virtual 
corrections contribute. We also included the box contribution (Fig. ||) into the "direct direct" part 
since it is known to be sizeable [Q, although it is formally a NNLO contribution. 

In the "direct fragmentation" part, the final state photon comes from the fragmentation of a 
hard parton participating in the short distance subprocess. From a technical point of view, a final 
state quark-photon collinear singularity appears in the calculation of the subprocess 717 — > jqq. At 
higher orders, final state multiple collinear singularities appear in any subprocess where a high pt 
parton (quark or gluon) undergoes a cascade of successive collinear splittings ending up with a quark- 
photon splitting. These singularities are factorized to all orders in a s and absorbed, at some arbitrary 
fragmentation scale Mp, into quark and gluon fragmentation functions to a photon, D 7 / C (z, M F ). 
When the fragmentation scale Mp, chosen of the order of the hard scale of the subprocess, is large 
compared to any typical hadronic scale ~ 1 GeV, these functions behave roughly as a/a s (M F ). 
Then a power counting argument tells that these fragmentation contributions are asymptotically of 
the same order in a s as the Born term. A consistent NLO calculation thus requires the inclusion of 
the 0{a s ) corrections to these contributions. 

Note that the singularity appearing in the process 717 — * jqq when the final state photon is emitted 
by the quark and becomes collinear, is subtracted and absorbed by the fragmentation function at 
the scale Mp, as explained above. Therefore both the "direct direct" and the "direct fragmentation" 
parts separately depend strongly on Mp and the attribution of the finite terms to either of these 
parts is scheme dependent. Only in the sum of these parts the Mp dependence flattens as expected. 

The collinear singularities appearing at NLO if the incident photon splits into a collinear qq 
pair are absorbed into the functions F|j/ 7 (a; 7 ,M) at the factorization scale M. (Analogous for the 

3 By "prompt" we mean that the photon is not produced from the decay of light mesons. 
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Figure 1: Examples of direct direct and direct fragmentation contributions at leading order 




rms~s~s~is^ 



Figure 3: The box contribution 



proton distribution functions F b / p (x p , M p ) ; we will set M p = M in the following.) Thus, by the 
same reasoning as above for the final state, the "initial direct" and "initial resolved" parts separately 
show a strong dependence on M which cancels out in the sum. Therefore it has to be stressed that 
only the sum over all four parts has a physical meaning. Figure [5] illustrates these cancellation 
mechanisms. 

The overall reduction of the scale dependence when going from leading to next-to-leading order 
can be seen in Fig. [|. The scales Mp and M have been set equal to /1, and /1 has been varied between 
[i = pip/2 and /1 = 2 pip. One can see that the NLO cross section is much more stable against scale 
variations, it varies by less than 10% in this /i range. 
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Figure 4: Dependence of the total cross section on scale variations, /i = M — Mp is varied between 
fi = Pp/2 and fi — 2p1p. 
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M F dependence, initial direct, final dir/frag 
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Figure 5: Cancellation of the leading dependence on the fragmentation scale Mf between con- 
tributions from direct and fragmentation final states, and on the factorization scale M between 
parts with direct and resolved initial state. The results are normalized to the total cross section at 
M F = M = [i=p 1 T . 
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2.2 Treatment of infrared singularities 

There are basically two methods to isolate the infrared singularities appearing in the calculation at 
NLO: The phase space slicing method and the subtraction method Jjl|. The method used here 
follows the approach of |ll], [l4) which combines these two techniques. We will outline the strategy 
only shortly, for more details we refer to |jjfl . 

For a generic reaction 1 + 2^3 + 4 + 5, at least two particles of the final state, say 3 and 4, 
have a high px and are well separated in phase space, while the last one, say 5, can be soft, and/or 
collinear to either of the four others. In order to extract these singularities, the phase space is cut 
into two regions: 

- part I where the norm px5 of the transverse momentum of particle 5 is required to be less 
than some arbitrary value prm taken to be small compared to the other transverse momenta. 
This cylinder contains the infrared and the initial state collinear singularities. It also contains 
a small fraction of the final state collinear singularities. 

- parts II a(b) where the transverse momentum vector of the particle 5 is required to have a 
norm larger than pTm, and to belong to a cone C^C^) about the direction of particle 3(4), 
defined by (775 — iji) 2 + (</>s — (pi) 2 < R 2 h (i = 3,4), with R t h some small arbitrary number. 
63(64) contains the final state collinear singularities appearing when 5 is collinear to 3(4). 

- part II c where PT5 is required to have a norm larger than PTm, and to belong to neither of 
the two cones C3, C4. This slice yields no divergence, and can thus be treated directly in 4 
dimensions. 

The contributions from regions I and IIa,b are calculated analytically in d = 4 — 2e dimensions and 
then combined with the corresponding virtual corrections such that the infrared singularities cancel, 
except for the initial (resp. final) state collinear singularities, which are factorized and absorbed into 
the parton distribution (resp. fragmentation) functions. 

After the cancellation, the finite remainders of the soft and collinear contributions in parts I 
and IIa,b,c separately depend on large logarithms lnpr m , In 2 pTm and \riRth- When combining the 
different parts, the following cancellations of the prm and Rth dependences occur: 
In part I, the finite terms are approximated by collecting all the terms depending logarithmically on 
PTm and neglecting the terms proportional to powers olpTm- On the contrary, the R t h dependence 
in the conical parts II a and II b, is kept exactly. This means that an exact cancellation of the 
dependence on the unphysical parameter R t h between part lie and parts II a, b occurs, whereas the 
cancellation of the unphysical parameter prm between parts II c, II a,b and part I is only approximate. 
The parameter prm must be chosen small enough with respect to p T in order that the neglected 
terms can be safely dropped out. On the other hand, it cannot be chosen too small since otherwise 
numerical instabilities occur. We have investigated the stability of the cross section by varying 
PTm and Rth between 0.005 and 0.1 (see Figure [|) and accordingly chosen the optimal values 
PTm = 0.05 GeV, Rth = 0.05. 
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Figure 6: Dependence of the total cross section on variations of the slicing parameter pTm- Ra 
denotes the total (nonisolated) cross section normalized to the total cross section evaluated with 
PTm — 0.005, R<j{pTm) = cr tot (pTm) / 'c toi (prm = 0.005). One can see that there is a plateau where 
the cross section is fairly insensitive to variations of pTm- The same study has been made for the 
dependence on R t h, but there the cross section is completely stable within the numerical errors since 
the R t h dependence has been kept exactly in all parts of the matrix element. 
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2.3 Implementation of isolation cuts 

In order to single out the prompt photon events from the huge background of secondary photons 
produced by the decays of 7r°,?y,w mesons, isolation cuts have to be imposed on the photon signals 
in the experiment. A commonly used isolation criterion is the following^: A photon is isolated if, 
inside a cone centered around the photon direction in the rapidity and azimuthal angle plane, the 
amount of hadronic transverse energy E T ad deposited is smaller than some value Ex max fixed by 
the experiment: 



(r, - Vl ) 2 + (0 - 7 ) 2 < i4 p 

£j t S max 



(4) 



Following the conventions of the ZEUS collaboration, we used ETmax = £Pr with e = 0.1 and 
-Rexp = !■ Isolation not only reduces the background from secondary photons, but also substantially 



reduces the fragmentation components, as will be illustrated in section 3.1. 

Furthermore, it is important to note that the isolation parameters must be carefully fixed in 
order to allow a comparison between data and perturbative QCD calculations. Indeed a part of the 
hadronic energy measured in the cone may come from the underlying event; therefore even the direct 
contribution can be cut by the isolation condition if the latter is too stringent. Let us estimate the 
importance of this effect and assume that the underlying event one-particle inclusive distribution is 
given by 

dn^ h 4 2p T 

= e <PT ' , (5) 

Pt dpT dri d(j) 2tt (px) 2 

being normalized to n particles per unit of rapidity. The probability that the isolation condition 
is fulfilled by a particle from an underlying event is 

« f AAA r A dn(1) 

n isoi = / dcpdrj prdpT -j— tt 

Jcone JE Tmax VT dpT df] d(f> 

With the ZEUS isolation parameters, ETmax — 0.5 GeV for a photon of p T = 5GeV. Using n — 3 
and (pt) ~ 0.35 GeV extracted from pS, one obtains 

« 0.33 

This estimation is very rough and underestimates the true effect because there is also a non- negligible 
probability to fulfill the isolation condition with two underlying particles falling into the cone. Only 
a detailed Monte Carlo description of the underlying events can allow a reliable estimate of this 
non-perturbative effect. Here we just note that the cut put by ZEUS (ETmax ~ 0.5 GeV) is likely 
to be too low to eliminate any underlying event contamination and therefore makes a comparison 
between the partonic level QCD predictions and the (hadronic level) data difficult. 

2.4 Features of the code 

The code consists of four subparts corresponding to each of the four categories of subprocesses. 



For each category, the functions corresponding to the parts I, IIa,b,c described in section 2.2 are 
integrated separately with the numerical integration package BASES (2^]. Based on the grid produced 
by this integration, partonic events are generated with SPRING |23| and stored into an NTUPLE or 
histogrammed directly. It has to be emphasized that we generate final state partonic configurations. 

4 A more sophisticated criterion has been proposed in [ |17[ , in which the veto on accompanying hadronic transverse 
energy is the more severe, the closer the corresponding hadron to the photon direction. It has been designed to make 
the "fragmentation" contribution vanish in an infrared safe way. 
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Hence this type of program does not provide an exclusive portrait of final states as given by hadronic 
event generators like PYTHIA or HERW1G. On the other hand, the latter are only of some improved 
leading logarithmic accuracy. The information stored in the NTUPLE are the 4-momenta of the 
outgoing particles, their types (i.e. quark, gluon or photon), the energy of the incident photon and, 
in the fragmentation cases, the longitudinal fragmentation variable associated with the photon from 
fragmentation. Furthermore a label is stored that allows to identify the origin of the event, e.g. if 
it comes from a2^2ora2->3 process. Based on the information contained in these NTUPLES, 
suitable observables can be defined and different jet algorithms can be studied. The isolation cuts 
are included already at the integration level, but the user of the program can turn isolation on or 
off and vary the input parameters for the isolation cut at will. 



3 Numerical results and comparison to ZEUS data 

In this section we present some numerical results for isolated prompt photon production. We restrict 
ourselves to the inclusive case, photon + jet production will be discussed in detail in a forthcoming 
publication. 

For the parton distributions in the proton we take the MRST2 |l9| parametrization. Our de- 
fault choice for the photon distribution functions is AFG |2(J, for comparisons we also used the 
GRV (2^] distributions transformed to the MS scheme. For the fragmentation functions we use the 
parametrization of Bourhis et al p2[| . We take nf = 4 flavors, the contributions from the b quark can 
be neglected since the typical energy scale of the partonic process is too close to the 6-production 
threshold. For ct s (p) we use an exact solution of the two- loop renormalization group equation, and 
not an expansion in log(/i/A). Unless stated otherwise, the scale choices M = Mf = fJ> = p\ have 
been used. The rapidities refer to the ep laboratory frame, with the HERA convention that the 
proton is moving towards positive rapidity. HERA operates with an electron energy of 27.5 GeV 
and a proton energy of 820 GeV. 

Following the ZEUS collaboration |D| , the parameters for the Weizsacker- Williams spectrum are 
Qmax = i GeV and the photon energy y = E 1 / E e has been restricted to the range 0.2 < y < 0.9. Note 
that experimentally, the energy of the incoming photon in photoproduction processes is reconstructed 
from the final hadron energies with the Jacquet-Blondel method, 

where the sum is over all calorimeter cells, E is the energy deposited in the cell and p z = E cos 0. In 
order to obtain the "true" photon energy y, corrections for detector effects and energy calibration 
have to be applied to yjs- These corrections are assumed to be uniform over the whole y range and 
enter into the experimental systematic error. However, as the background varies with the photon 
energy y, these corrections may not be uniform. It has to be emphasized that the cross section is 
very sensitive to a variation of the energy range of the photon. (See Figure [L2 and discussion below.) 



3.1 Numerical results for inclusive prompt photon production 

If not stated otherwise, all plots showing the photon rapidity {rf 1 ) dependence are integrated over 
5 GeV < pi < 10 GeV and 0.2 < y = E y /E e < 0.9. 

Figure H shows a comparison of the NLO to the leading order result for the isolated cross section 
da I 'dtp The importance of the box contribution is clearly visible. The higher order corrections 
enhance the isolated cross section by about 40%. 

Fig. U shows the rapidity distribution of the full cross section before and after isolation. As 
already mentioned in section 2.3, we used the isolation cuts Exmax = £Pr with e = 0.1 and R, 



cxp 



1 to match those of the ZEUS collaboration. Fig. H also shows the effect of isolation on the 
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fragmentation partj^] separately. Isolation reduces the fragmentation component to about 6% of the 
total isolated cross section. 

In Fig. ^ the relative magnitude of all four components contributing to da ep ^ lX / 'drp before and 
after isolation is shown. Note that isolation increases the contributions with a direct photon in the 
final state slightly since there the cut mainly acts on a negative term, which is the one where parton 
5 is collinear to the photon. It should be emphasized that Figure |9] has to be read with care since 
the individual parts have no physical meaning and are very sensitive to scale changes. Nevertheless 
the dominance of the resolved direct part remains if we choose e.g. /j, = M = Mp = pJ/2 or 2 pj. 

Figure [l^ shows the relative magnitude of contributions from resolved and direct photons in 
the initial state to the isolated cross section. From the pj, distribution one can conclude that the 
resolved part dominates the cross section for small values of jP T such that it would be useful to look 
at the photon rapidity distribution at = 5 GeV in order to discriminate between different parton 
distribution functions in the photon. 

Since the gluon distribution in the photon is of particular interest, the sensitivity to the gluon in 
the photon is investigated in Fig. [n]. One can see that the gluon distribution in the photon starts 
to become sizeable only for photon rapidities rp > 1 and dominates over the quark distribution for 
about rp > 2.5 . Therefore the region of large photon rapidities and small photon px is the one where 
the sensitivity to the gluon in the photon is largest. In order to test further the sensitivity to the 
gluon, we increased the gluon distribution in the photon uniformly by 20%. As can be anticipated 
from Fig. (ill the effect becomes sizeable only for n 1 > 2 and leads to an increase of the cross section 
by about 10% for rf > 2.5 . We conclude that in the region rf 1 < 1, there is basically no sensitivity 
to the gluon in the photon. However, investigating the photon + jet cross section instead of the 
inclusive case offers larger possibilities to constrain the gluon in the photon since there one can vary 
the photon and the jet rapidities in order to single out a kinematic region where the sensitivity is 
large ||. 

Figure 12 shows the effect of a ten percent uncertainty in the "true" bounds of the photon energy 
y. One can see that a change of the lower bound on y has a large effect, in particular at large photon 
rapidities. This comes from the fact that the Weizsacker- Williams distribution is large and steeply 
falling at small y. Increasing the lower bound on y therefore removes a large fraction of the direct 
events with lower energy initial photons, (y = x e for the direct events and large rp correspond 
to small x e .) At large photon rapidities the spread due to the use of different parton distribution 
functions for the photon is smaller than the one caused by a 10% variation of the lower bound on 
y. On the other hand, the region of large photon rapidities is of special interest since there the 
gluon in the photon is dominating. Therefore a small experimental error in the reconstruction of the 
"true" photon energy is crucial in order to be able to discriminate between different sets of parton 
distribution functions in the photon. 

It has been tested that the effect of using different proton distribution functions - for example 
the CTEQ4M or the MRST1 set of proton distribution functions - is of the order of 3% at most. 
In all photon rapidity bins the spread is smaller than the one caused by different sets of photon 
distribution functions (which is about 10% at small photon rapidities, see e.g. Fig. Jbi]). Thus a 
discrimination between different sets of photon distribution functions should be possible with the 
forthcoming full 1996-2000 data set analysis, where the errors on the data are expected to be small 
enough. 



5 By "fragmentation part" we mean the sum of the two subparts "direct fragmentation" and "resolved fragmenta- 
tion" . Analogously, " resolved" denotes the sum of " resolved direct" and " resolved fragmentation" . 
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photon rapidity and transverse momentum distribution, with isolation. 
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Figure 12: Photon rapidity distribution da ev ^ lX j 'drf for isolated prompt photons integrated over 
5 GeV < pip < 10 GeV and different lower bounds on y. Solid line: 0.2 < y < 0.9 with AFG photon 
structure functions, dotted line: bounds on y changed by about 10%, dashed line: 0.2 < y < 0.9 
with GRV photon structure functions 
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3.2 Comparison with ZEUS data 



In this section we compare our results to the ZEUS 1996-97 data on inclusive prompt photon 
photoproduction |ll| . Figures [l3| and |lj show the photon pt and rapidity distributions with AFG 
resp. GRV sets of structure functions for the photon. For the px distribution the agreement between 
data and theory is quite good. In the rapidity distribution (Fig. |lj) the data fluctuate a lot, such 
that the agreement is still satisfactory. However, it seems that theory underpredicts the data in the 
backward region, whereas the theoretical prediction tends to be higher at large photon rapidities. 
The curves of Gordon and Krawczyk/Zembrzuski || given in also show this trend. At high 
rf 1 the reason for the difference could be that the isolation cut in the experiment removes more 
events than in the theoretical (parton level) simulation, as discussed in section 2.3. 

Figure |l5| shows that the discrepancy between theory and data at low rf 1 comes mainly from the 
domain of small photon energies, whereas the discrepancy at large rf is only present in the range 
of large photon energies. Note that at large rp and large y the resolved part dominates and the 
underlying event could have a large multiplicity. Therefore the isolation criterion could also cut on 
the non- fragmentation contributions as discussed in section 2.3. 
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Figure 13: Comparison to ZEUS data of photon px distribution da ep ^" <x / dp T for isolated prompt 
photons. Results for two different sets of parton distributions in the photon are shown. 
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comparison to ZEUS data, with isolation 
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Figure 14: Comparison to ZEUS data of photon rapidity distribution d<r ep ^ lX jdrp for isolated 
prompt photons. 
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subdivisions of photon energies 
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Figure 15: Photon rapidity distribution da ep ^ lX /drp integrated over 5 GeV< Pj. <10 GeV and 
different subdivisions of photon energies: (a) 0.2 < y < 0.32, (b) 0.32 < y < 0.5, (c) 0.5 < y < 0.9. 
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4 Conclusions 



We have presented a program for prompt photon photoproduction which includes the full next-to- 
leading order corrections to all contributing subparts. It is a general purpose code of partonic event 
generator type and as such very flexible. 

We used it to study the possibility to constrain the quark and gluon distributions in the photon. 
It turned out that the sensitivity to the gluon distribution in the photon is negligible in the rapidity 
range -0.7 < rp < 0.9 studied by ZEUS. 

A discrimination between the AFG/GRV sets of parton distributions in the photon is not possible 
with the present experimental errors on the ZEUS 1996/97 data. However, a forthcoming analysis 
of all 1996-2000 data announced by the ZEUS collaboration will drastically improve this situation. 

We have shown that the cross section is very sensitive to small variations of the photon en- 
ergy range. Therefore a good control of the experimental error on the photon energy fraction y 
(reconstructed experimentally from the Jacquet-Blondel variable ujb) will be crucial for future 
comparisons. 

Despite the large fluctuations of the data, one can say that there is a trend that theory overpre- 
dicts the data in the forward region. The reason might be that the isolation cut imposed at partonic 
level in the perturbative QCD calculation does not have the same effect as the experimental one. If 
the experimental cut is too stringent, a large fraction of the hadronic energy in the isolation cone 
may come from underlying events, such that experimentally a larger number of events is rejected. 
We gave a rough estimate of the underlying events to be expected in the isolation cone. 

The possibilities offered by the study of photon + jet photoproduction will be investigated in a 
forthcoming publication p4| . 
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